Griffiths singularity in the random Ising ferromagnet 
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in ■ 

■ The explicit form of the Griffiths singularity in the random ferromagnetic Ising model in external 

' magnetic field is derived. In terms of the continuous random temperature Ginzburg-Landau Hamil- 

£nQ I tonian it is shown that in the paramagnetic phase away from the critical point the free energy as 

the function of the external magnetic field h in the limit h — > has the essential singularity of the 
form exp | — (const) /h ^ 3 ^ (where 1 < D < 4 is the space dimensionality). It is demonstrated that 

in terms of the replica formalism this contribution to the free energy comes due to off-perturbative 
replica instanton excitations. 

I. INTRODUCTION 

The history of the problem of the Griffiths singularities starts from the theorem of Lee and Yang which states 
that the partition function of an (ordered) ferromagnetic Ising model (in any space dimensions and for any lattice 
connectivity) as the function of the external magnetic field h in the thermodynamic limit has a continuous distribution 
of zeros along the imaginary axis of the complex parameter h — x+iy. Moreover, in the paramagnetic phase depending 
on the temperature this distribution starts at finite distance Ay from the real axis, which means that here the free 
energy of the system is an analytic function of the real magnetic field h (with y — 0). On the other hand, when the 
temperature T approaches the phase transition point T c from above the value of the interval Ay shrinks to zero, so 
, that the distribution of zeros touches the real axes at T = T c , which indicates (in agreement with the modern theory 
£5 ' of the critical phenomena) that at the phase transition point the free energy must be a non-analytic function of (real) 
ft-»0. 

i Now, if one would try to apply the above general observations for random systems, the consequences turn out to 

be much more tricky. As an example, let us consider bond diluted ferromagnetic Ising model, in which the critical 

t-H , temperature T c (p) is the (decreasing) function of the degree of dilution p. Here one can note that in an infinite 
system at a temperature T which is above T c (p) (such that the state of the system is paramagnetic) but below the 
critical temperature T of the corresponding pure system, with a finite (exponentially small) probability there are 
exist arbitrary large less diluted ferromagnetic "islands" which are critical exactly at this given temperature T. It 

ly-^ \ is important that such clusters exist at any temperature in the interval T c (p) < T < Tq. Thus one can expect that 
. the free energy of such random system must be a non-analytic function of the external magnetic field ft, — * at any 
' temperature between T c (p) and To, Ref 2]. Unlike pure systems, however, here it is much more difficult to predict the 

<^> explicit form of such non-analyticity. 

For the one-dimensional diluted Ising chain it has been shown that its free energy F(T, h) is non-analytic in the 
point T — 0, h = 0, and moreover the divergences of the pure system thermodynamics are replaced by an essential 
singularity at which all functions are finite and infinitely differentiable- . According to further more general studies 
1 , \ (in terms of heuristic arguments** and the Bethe lattice systems*) , the form of this singularity has been argued be 
of the type exp [— (const) /h]. 

At the same time there has been much interest in the dynamical properties of such systems in the temperature 
interval T c (p) < T < Tq. It has been discovered that the relaxation processes here are slower than the exponential 6 . 
It turned out that due to the presence of rare large ferromagnetic clusters the relaxation of e.g. the order parameter 
takes either " stretched-exponential" form exp [— (t/r) l3 ^ 1 \ j controlled by the temperature dependent exponent j3 < 
1 (the result preferred by the numerical simulations* ), or even slower type of decay exp [— (const) (int) ^ ^ 1 ^ 
(predicted analytically 6 **) . Due to these quite non-trivial dynamical properties, which are essentially different from 
the paramagnetic phase, the state of the system in this temperature interval is usually called the Griffiths phase. 

Besides dynamics, an essential progress has been achieved in the analytical investigation of the distribution of 
zeros of the partition function alone the imaginary axis of the complex magnetic field. In particular, formal replica 
supersymmctric calculations performed for the random temperature Ginzburg-Landau Hamiltonian 9 has demonstrated 
the importance of the instanton- like off-perturbative excitations which provide the development of the "tail" in the 
distribution of zeros in analogy with the density of states in the Anderson localization problem. On the other hand, 
the study of the diluted infinite-range Ising ferromagnet with finite connectivity in purely imaginary magnetic field 
h = iy has shown that in the paramagnetic phase the tail of the density of zeros p(y) has the explicit form of the type 
p(y) ~ exp [— f(T)/y], where the function of the temperature f(T) vanishes in the critical point 10 . Similar tail has 
been also found numerically in the two-dimensional diluted Ising model*k All that, however, doesn't answer the main 
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question, what is the explicit form of the Griffiths singularities of the thermodynamical quantities as the function of 
the real magnetic field in the point h — 0. 

In this paper I am going to consider this problem in terms of the Z?-dimensional random temperature Ginzburg- 
Landau Hamiltonian (in dimensions 1 < D < 4) in the paramagnetic phase away from the critical point. In the next 
section simple heuristic arguments will be proposed which demonstrate on a qualitative level the physical mechanism 
by which non- analytic contributions appear in the thermodynamical functions of such type of systems. After that, the 
systematic method of off-perturbative replica calculations will be formulated in Section III. Finally, in Section IV it 
will be demonstrated that non-analytic (Griffiths) contributions to the thermodynamics comes from off-perturbative 
instanton-like excitations. In the limit h — > such contributions to the free energy is argued to have the explicit form 



AF 



exp 



-C h- D ' z 



(1.1) 



where the constant C is defined by the parameters of the Ginzburg-Landau Hamiltonian and by the strength of the 
disorder. 



II. HEURISTIC ARGUMENTS 



Before starting doing systematic calculations let us try to understand on a pure qualitative level, using simple 
" hand-waving-arguments" , what is the physical mechanism by which non-analytic contributions are coming into the 
free energy (and others thermodynamical functions) in the paramagnetic phase of weakly disordered ferromagnetic 
Ising model in external magnetic field. Let us suppose that such system in D dimensions can be described by the 
continuous Ginzburg-Landau Hamiltonian 



H = d u x 



l - (V0(z)) 2 + \ (r - Sr(x)) 2 (x) + +ig0 4 (x) - h<f>(x) 



(2.2) 



where the reduced temperature parameter r is taken to be positive and sufficiently large to place the system into 
the paramagnetic phase. The disorder is modeled by a random function t(x) which is described by the Gaussian 
distribution, 

P[St] = p expf-^ J d D x{8r{x))^j , (2.3) 

where u is the small parameter which describes the strength of the disorder, and po is the normalization constant. 

Intuitively, it is clear that non-trivial contributions to the thermodynamics are coming due to rare "ferromagnetic 
islands" in which the value of St(x) is on average bigger than r. Let us consider such an island, which is characterized 
by the linear size L and the typical value of the "local temperature" (r — St) = — £ < Its probability is exponentially 
small, 

and therefore such islands are well separated from each other and can be considered as non-interacting. Note that the 
physical machanism of slowing down of the relaxational processes due to the presence of these ferromagnetic islands 
are more or less clear. The magnetization orientation of the ferromagnetic cluster can be either "up" or "down", and 
if its size L is big, then it would require a long "elementary relaxation time" t(L) to flip it from one orientation to the 
other (since flipping of the cluster would require overcoming a big energy barrier, which is proportional to the volume 
of the cluster L D ). The origine of the non-analytic contributions appearing in terms of pure statistical mechanics is 
much less clear: since time is formally infinite here, the presence of big energy barriers, separating the two orientations 
is irrelevant. 

First, let us consider what is going on in the zero external magnetic field. Here one could distinguish two types of 
contributions to the thermodynamics: 

(1) the perturbative one, coming from the spatial scales smaller than the correlation length R c , which formally 
could be computed e.g. in terms of the renormalization-group (RG) approach-^; 
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(2) the non-perturbative contributions (missing in the RG treatment) due to the "up" and "down" ferromag- 
netic states of rare ferromagnetic islands discussed above, which are coming from the spatial scales bigger than the 
correlation length R c . 

It has to be noted that the island with small (negative) local temperature — £ can be characterized as having 
the distinct (mean-field) "up" and "down" states only if its size is much bigger than its local correlation length 
Rc(0 ~ i~ x ^ 2 In what follows we are going to use the mean-field values of the critical exponents assuming that 
the temperature of the system is taken sufficiently far away from the critical region near T c (the size of this region 
is of the order of r g — g 2 /^~ D > which is small if the coupling parameter g is taken to be small). Thus, with the 
exponential accuracy the (Griffith-like) contribution to the free energy coming from these local ferromagnetic states 
can be estimated by their probability, ea. (|2.4|l : 



G 



oo poo 

d£ / dL exp 



dt; exp 



(const) 



-D/2 



Here in the integration over £ the leading contribution comes from the vicinity of the saddle-point value 

D 



4- D 



(2.5) 



(2.6) 



(which is positive in dimensions D < 4, and r g provided r 3> r g ). In this way we obtain the following estimate 
for the off-perturbative contributions coming from rare locally ferromagnetic islands: 



Fq ~ exp 



— (const) - 



r (4-D)/2 



(2.7) 



In fact, this result (including the value of the (const) factor), as we will see in section IV, can be derived analytically 
in terms of the formal replica calculations as the contribution from the localized (instanton-like) solutions of the 
mean-field saddle-point equations^. 

In the presence of non-zero external magnetic field h the situation becomes slightly more tricky. Let us consider 
again the ferromagnetic island of the size L with negative "local temperature" (r — St) = — £ which is described by 
the Hamiltonian, ea. H2.2|) . Note first of all, that the effective potential 



has two minima only if the value of £ is not too small, namely at 



(2.8) 



(2.9) 



otherwise, at £ < it has a unique minimum (see Fig.l). In the case of the two minima, the state "up" (4> > 0) 

has the lower energy E u (£, h, L), and it is the ground state of the island, while the state "down" (<f> < 0) can be 
considered as the excitation since it has higher energy, Ed(t;,h, L). The crucial point is that the excited state, Ed 
disappears discontinuously at £ = (h) such that the energy difference AE = Ed — E u remains finite at this critical 
value of £. Thus, the contribution to the total free energy of such ferromagnetic islands can be estimated as follows: 



F G (h) 



dt 

t.(h) JR c (0 



dLV[L,£] in [e- E - + e- Ed ] 



d£ / dL V [L, £] ln[e _B »] (2.10) 
>Rc(t) 



where V[L,£] is the probability to have the island of the size L with the local temperature £, ea. H2.4|) . Since the 
integration over L sticks to the lower bound R c , we find: 



F G (h) 



d£P(t;)E u (Z,h) - / d^P(0 In 



l + e 



-AE(£,h) 



where 



P(0 ~ exp 



-(const) 



(r + OV 



D/2 



(2.11) 



(2.12) 



(a) (b) (c) 

FIG. 1: Qualitative shape of the potential 17(0), eq.JUHJ): (a) at £ > £,(/i); (b) at £ > £*(h); (c) at £ = £*(h). 

On the other hand, similar considerations for the zero field case (when Ed = E u ) yield 



F G (0) 



OO POC 

d£ / dLV[L,£] In 



2e -B„(e,o) 



<%P(0 [£ u (£,0)-(ln2)] 



(2.13) 



Thus, for the free energy difference, AFg(/i) = Fo(h) — Fg(0) we find: 



AF(h) 



<%P(Z) [E u (Z,h)-E u (£,0)] 



(h) 



d£ P(0 In 



1 + e 



(ln2) / d£ P(0 (2.14) 



In the limit /i — > the first two terms in the above equation provide regular functions of h (both (_£?„(£, h) — E u (£, 0)) 
and AE(£, h) go to zero as a power functions of h in the limit h — > 0), while the last one is just the non-trivial Griffiths 
contribution 5Fq(K) which has the form of the essential singularity: 



SF G (h) - (ln2) J d£ P(0 - exp (-(const)^* D/2 (h)\ 
Substituting here the value £*(h) — h 2 /' i g 1 / 3 1 ea. (|2.9fl . one eventually finds: 



(2.15) 



SFcih) ~ exp —(const) 



3 D/6 



h-°/ 3 



(2.16) 



In Section IV it will be demonstrated how this result can be derived in terms of the formal replica calculations. But 
first we have to formulate the general lines of the replica approach for the off-perturbative contributions. 



III. OFF-PERTURBATIVE REPLICA CALCULATIONS 



In this section I am going to formulate a general systematic approach for the calculations of off-perturbative 
contributions (if any) coming from local minima states, which in the configurational space are well separated from 
the ground state. 

Let us consider a general random system described by a Hamiltonian H [<p(x)], and let us suppose that in addition 
to the ground state, there is another thermodynamically relevant (Griffith) region of the configurational space located 
"far away" from the ground state and separated from it by a finite barrier of the free energy (see Fig. 2). In other 
words, we suppose that the partition function (of a given sample) can be represented in the form of two separate 
contributions: 



Z = I V4>{x)e-P H = e-P F ° + e-^ = Z Q 



(3.17) 
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FIG. 2: Schematic structure of the free energy landscape in the case of two well separated thermodynamically relevant valleys 
of the configurational space 

where Fq is the contribution coming from the vicinity of the ground state, and F± is the contribution of the Griffiths 
region. Then, for the averaged over disorder total free energy we find: 

1, — 1: 



T=--\nZ = F Q --\n[\ + Z 1 Z^ 1 ] (3.18) 
The second term here, which is just the Griffiths contribution, can be represented as follows: 

F n =-- V zrZ- m = ~ lim V Z n (m) (3.19) 

8 ^— ' to ' ft n^o ^ to 

m—l 



where 



z n (m) = n / n / e-^-[*? ) --^ ) ^° , .-.e ra ] ( 3. 20) 

6=1 ^ c=l ^ 

is the replica partition function (H n [<p] is the corresponding replica Hamiltonian), in which the replica symmetry in 
the n-component vector field <p a (a — 1, ...,n) is assumed to be broken. Namely, it is supposed that the saddle-point 
equations 

5H n [4>] 

— — — = 0, (o=l,...,n) (3.21) 

0(pa{x) 

have non-trivial solutions with the RSB structure 

{4>i(x) for a = 1, to 
(3.22) 
4>o(x) for a = m + 1, ...,n 

with 0i (x) 7^ 00(^)7 so that the integration in the above partition function, ea. (|3.2U() . goes over fluctuations in the 
vicinity of these solutions: 

4>f\x) = (f>i(x) + ifb(x), (b=l,...,m) 

4>W(x) = Mx)+Xc(x), (c=l,...,n-m) (3.23) 

It should be stressed that to be thermodynamically relevant, the RSB saddle-point solutions, ea. (|3.22|l . should 
satisfy the following tree crucial conditions: 
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(1) the solutions should be local in space, so that they are characterized by finite space sizes R(m); in this case 
the partition function, ea. l|3.2U|l . will be proportional to the entropy factor V/R D (m) (where V is the volume of the 
system), and the corresponding free energy contribution Fq, ea. H3.19(l . will be extensive quantity; 

(2) they should have finite energies E(m) = H n [</>*]; 

(3) the corresponding Hessian matrix of these solutions should have all eigenvalues positive. 

Thus, in the systematic calculations one should find all saddle-point RSB solutions <j>a(x) (satisfying the above three 
requirements), ea. l|3.22|l . after that one has to compute their energies E(m) (for n — > 0), next one has to integrate 
over the fluctuations in the vicinity of these solutions, and finally one has to sum up the series 



V 



(-1) 



m — 1 



-i? 



— D 



I ~\-l/ 2 

(to) detT 

V / n=0 



(3.24) 



where T is the (n x n) matrix 



T i — 

± an/ — 



S 2 H [4 



(3.25) 



The above scheme of calculations can be easily generalized for an arbitrary number of the Griffiths regions. For 
example, let us consider the situation which is qualitatively represented in Fig. 3, when in addition to the ground state, 
the system has two thermodynamically relevant Griffiths states (which is just the case for the considered random Ising 
model) . In this case instead of eq. (|3.17|l we will have 




configurational 
space 



FIG. 3: Schematic structure of the free energy landscape in the case of three well separated thermodynamically relevant valleys 
of the configurational space 



Z = / V(j)(x)e 



-I3H _ P -I3F Q , P -$F X , P -/3F 2 



an! correspondingly, instead of ea. l|3.19|l we find 

1- 



= Z + Zi 



Fa 



{-iy 



i v 

— lim > 

3 n^Q *-~> 



(-1) 



m-1 m 



m— 1 



fc=0 



Zo 



(3.26) 



J2c? (z?z™~ k z- m ) 



k=0 



(3.27) 
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where C™ = m!/(fe!(m — k)\) is the combinatoric factor. Here, in the replica partition function 

k „ rn—k n—m /> 

z n (k,m) = n / ^ n / n / ^ (o) e-^K 1 ^.^] 

the integration is supposed to be performed in the vicinity of the saddle-point replica vector 

4>i(x), for a = 1, k 

4>* a { x ) — { 02 (a?) > for a = fc + 1, to 
</>o (x), for a = m + 1, ...,n 



(3.28) 



(3.29) 



(where 4>i(x) ^ faix) ^ 4>o{x)) which is the solution of the saddle-point equations l|3.21ll . Finally, for the Griffiths 
contribution, instead of ea. H3.24[l one obtain 



Fr, 



(-1) 



m— 1 



m— 1 



0m 



J2 C T R~ D (detTj _ Q e 



-j3E(k,m) 



(3.30) 



fc=0 



where i?(fc, m) = H n — i-o [0*] is the energy of a given solution, ea. (|^.29|l . and T is the Hessian matrix, ea. H3.25(l . 

It is worthing to note that one can arrive to the same representations for the off-perturbative free energy contribu- 
tions, eas. (|3.24|) and (|3.3UI) . in terms of the so called vector replica symmetry breaking schem o 14 ! 15 , starting from the 
standard replica approach for random systems (F — — \im n ^o(Z n ~ l)/ n )- 

In the next section we will implement the programme described above for the concrete case of weakly disordered 
ferromagnetic Ising model in the high temperature paramagnetic phase. 

IV. REPLICA INSTANTONS 

After the standard Gaussian averaging over random 5t(x) (described by the distribution, ea. (|2.3|0 of the n-th power 
of the partition function one obtains the following replica Hamiltonian 



TL Tt TL Tt Tt 

L a— 1 a— 1 a— 1 a,b— 1 a— 1 

The corresponding saddle-point equations are 

-A<j> a (x) + T<pa(x)+g4>l(x)-ucf> a (x) I J2<t>b(x)} = h 



sf)=l 



Substituting here the anzats, ea. (|3.29(l . the above saddle-point equations are reduced to 

-A(j>i + T<j)i + g§\ — ucj)iS = h 

(i = 1, 2, 0) where 

n 

S = ^2<pl(x) = k(j)\ + {m - k)4> 2 2 + {n - m)4>l 



a=l 



which in the limit n — > turns into 



S = k<f>i + (to — k)4>\ — mcf>Q 



(4.31) 



(4.32) 



(4.33) 



(4.34) 



(4.35) 



Substituting eas. l|3.29|l . into the Hamiltonian, ea. l4.31|l . for the energy of this configuration (at n — » 0) we obtain: 



E{k, to) 



d D x 



2 ( V 0O + ( v 02) 2 



^(V0o) 2 + C/(0i,0 2 ,0o) 



(4.36) 
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where 



U((j>u<h,M = \ T [k<t>i + ( m - k )4>l - m<j)l] + jg [k<j)\ + (to - k)4>\ - m$j] 



u \k<t>\ + (to — k)<p2 — Trufaol — h \k<f>x + (m — k)(f>2 — to^o] 



(4.37) 



and the fields 4>i (x) , 4>2 {%) and <f>o(x) are defined by the equations (|4.33|l and 14.35fl . 



A. Zero external magnetic field 

First of all, we note that at h = 0, due to the symmetry (j) — > —<f> the solution of eas. (|4.33|l - 14.35|l takes the form: 



(j>x(x) = ~4>2(x) = 4>{x) 
<j, Q {x) = 

where the function <f>{x) is defined by the equation 

-A4>(x) + t4>(x) - \(m)(f>(x) 3 = 

which is controlled by the parameter 

A(m) = ura — g 



(4.38) 
(4.39) 
(4.40) 



In what follows this parameter will be assumed to be positive. In other words, the solution, which we are going to 
derived below, exists only for m such that 



m > 



(4.41) 



Substituting eas. H4.38|) into ea. (|4.3ti|) - (|4.37|) for the energy of this solution we obtain 



E{k,m) = E(m) = m / d u 



(4.42) 



Note here, that one should not be confused by the "wrong" sign of the </> 4 term in the above expression (which for 
the usual field theory would indicate its absolute instability). In fact, as we will see below, the integration over 
the replica fluctuations around considered solution in the limit n — > yields the Hessian matrix which has all the 
eigenvalues positive (this is quite standard situation for the replica theory: in the limit n — > 0, when the number of 
certain parameters become negative, everything turns " up down" , so that minima of the physical quantities turns 
into maxima of the corresponding replica quantitiesi^i^) . 
Rescaling the fields, 



A(m) 



iP(x/R c (t)) 



(where R c (t) — r 1 ^ 2 ), instead of eq. 



one get the differential equation which contains no parameters: 
-Ai/>(z)+ip(z) -ip 3 {z) = 
Correspondingly, for the energy, ea. H4.42H . one obtains: 



E{m) = 



um — g 



-<.*-D)/*Eo 



where 



E = / d D z 



(4.43) 



(4.44) 



(4.45) 



(4.46) 
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The equation 14.44fl is well know in the field theory (see e.g. Ref . [l7| ) : in dimensions 1 < D < 4 it has the smooth 
(with ip'(O) = 0) spherically symmetric instanton-like solution such that: 

ip(z < 1) ~ ip(0) ~ 1, 

ip(z > 1) - e~ z -> 0. (4.47) 

The energy Eg, ea. (|4.46|l . of this solution is a finite and positive number. In particular, at D = 3, i[)q ~ 4.34 and 
i? — 18.90. Note that according to the rescaling, ea. H4.43|l . the size of this instanton solution in terms of the original 
fields <j>(x) is R c = t" 1 / 2 . This size does not depends on k and to, and it coincides with the usual mean-field correlation 
length of the Ginsburg-Landau theory. Note also that due to the symmetry <f> — > — <j> of the considered solution, its 
parameters do not depend on k. Thus, we can perform the summation over k in the series, ea. (|3.30[) . which yields the 
trivial factor 2 m : 



F G ^ -VR c (t)- d Y ( 1)T " - 2 m (dett) expf-£ ^^T( 4 - g >/ 2 ) (4.48) 
^-r J , , m V /n=o \ um — g ) 



(—\\ m ~ l / ,\-l/2 

m>[g/ii] 

In other words, the considered two-step structure, ea. l|3.29ll - <|3.30[l . is equivalent to the one-step anzats, ea. (|3.22|) . 



JtHu)^^) for a =1,..., to 
■>»■ { V (4.49) 
for a = to + 1 , . . . , n 

which has additional degeneracy factor 2 m . 

The final step is the integration over fluctuations, which define the Hessian factor (detT). Introducing small 
fluctuations <p a (x) near the above instanton solution, <fi a (x) — 4>* a {x) + tp a {x), in the Gaussian approximation we get 
the following Hamiltonian for the fluctuating fields: 



H[(p] ~ J d D :< 



n n 

{x)ip a fa' 

a—1 a,a' — l 



(4.50) 



where the matrix K aa i (x) contains the m x to block: 

K$\x) = (l- ^-^-^{x^)\ S bb — ^ 2 (^) (4.51) 

\ um — g J um — g 

(b, b' = 1, to) and the diagonal elements for the remaining (n — to) replicas: 

K%- m) = ( 1 - </> 2 (xVT)) S cc> (4.52) 

\ um — g J 

(c, d = m+ 1, n). To obtain the explicit result for the Hessian, let us approximate the instanton solution, ea. H4.47|l . 
by the #-like function 

f ^o, for < z < 1 
il>(z) = { (4.53) 
[ 0, for z > 1 

Then, the approximate Hamiltonian for the fluctuating fields takes much more simple form 



H[<p] * ^ E / Fi^+^^tpWM + oE / T^j/l^WI 2 (4.54) 

tt,o'=l J \P>^ \ ' a=l ■'IPKV^ 

Here the p-independent matrix K aa > is given by eas. (|4. 5111 - 14. 5211 . where instead of the function ^{xy/r) one has to 
substitute the constant ip . 

The integration over the modes with momenta p <C \fr (corresponding to the scales much bigger than the size of 
the instanton), which are described by the second term of the above Hamiltonian, gives the contribution of the form 
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exp(— nVfus)- This contribution is irrelevant in the limit n — > 0, which is quite natural, because all such contributions 
must be already contained in the perturbative part of the free energy F , ea. l|3.18[) . which we do not consider here. 
The integration over the modes with momenta p 3> \pr is slightly cumbersome but straightforward: 



(det f\ 



-1/2 



exp 



-D/2 



d D p 



Tr In (p 2 5 aa ' + rK a 



2 J p>>VT (2*Y 

The matrix under the logarithm in the above equation contains (m — 1) eigenvalues: 

2 / um - 3g 2 
Ai = p + t [ 1 - — r tp 



um — g 



one eigenvalue: 



and (n — m) eigenvalues: 



A 2 = P 2 + t 1 - 



um — 3g 2 



um — g 



2um 
um — g 



A 3 = p l + r 1 



um — g 



which are all positive. Substituting these eigenvalues into ea. H4.55|) . after simple algebra in the limit n 
eventually obtains the following result: 



( det TJ ~exp 
Finally, substituting this value into ea. (|4.48fl we find 



3m , n 

-Ml 



2{um — g)' 



(4.55) 



(4.56) 



(4.57) 



(4.58) 



one 



(4.59) 



Fa 



-VRc(t) 



-D 



oo 

E 

m>[g/u] 



(-1)' 



■;?) 



2 m exp 



um — g 



r (4-Z»)/2 + 



3 m 



2(um — g) 



# 2 



Here one can note that under condition 



t>t„ 



.2/(4-o) 



(4.60) 



(4.61) 



the second term in the exponential of eq. I|4.60f) (which is the fluctuations contribution) can be neglected compared to 
the first one. This is not surprising because eq. (|4.61|) is nothing else, but the good old Ginzburg-Landau condition 
which defines the temperature region away from T c , where the critical fluctuations are irrelevant, and the behavior of 
the system is described by the mean-field exponents. 

The exact summation of the series in ea. (|4.60() is rather tricky problem, but with the exponential accuracy it can 
be estimated in a very simple way. One can easily see that in the limit of weak disorder, at u <C g, the leading 
contribution in this summation comes from the region m 3> g/u 3> 1 (where the exponential factor in ea. (|4.60() 
becomes m-independent) and this contribution is 



F G - exp -Ei 



T (4--D)/2 



(4.62) 



We see that this result nicely coincides with the naive " hand- waving" estimate, ea. (|2.7|) ■ where the (const) factor is 
the instanton energy Eq (in three dimensions Eq ~ 18.9). 



B. Non-zero external magnetic field 



Technically, the situation in non-zero magnetic field becomes much more cumbersome, but on a qualitative level, 
the main idea of the approach remains very simple. According to the physical discussion of Section II, non-analytic 



11 



contribution to the free energy in the presence of external magnetic field appears due to the fact, that some of the 
instanton-likc configurations (of the type, considered in the previous subsection) disappear via a finite jump. 

To understand that, let us consider first the structure of the potential energy J7(0i, 02, 0o), ea. (|4.37[l . The extrema 
of this potential are defined by the three equations 

Tfa + g$ - ufcS = h (4.63) 

(i = 1, 2, 0), where 

S = k(j>l + {m- k)4> 2 2 - m4>l (4.64) 
Simple algebraic manipulations reduce these equations to 

0i + 02 + 0o = (4.65) 



and 



r 



0i<M0i+0 2 ) = -- (4-66) 
9 



+ [g + u(m - fc)] 4>{ + (g + uk) 02 + {g + 2um) 0i0 2 = (4.67) 



First, treating the magnetic field h here as a small correction to the zero-field solution, 



Ah=0) 
y 



(4.68) 



one easily finds 



gcp^ym) 
0(m) + O (h 2 ) , 

h 



g<p 2 (m) 



O (h 2 ) (4.69) 



It is clear that this linear in h shift of the extremum of the potential U (0i, 02, 03) provide not more than a linear in h 
corrections to the zero field instanton space configuration as well as to its free energy contribution considered in the 
previous subsection. It has to be noted, however, that the above result, ea. (|4.69() . is valid only until the summation 
parameters k and m are not too large: 

r /o\ 2 / 3 

m > k < - (t ( 4 - 7 °) 

Simple analysis shows that as k and m grow, the values of (negative) 02 (fc, m) and 0o(fc, m) become closer and closer to 

each other, and finally, one arrive to the critical cc 
into equations (|4.65|) - (|4.67|) one easily finds that 



each other, and finally, one arrive to the critical configuration when their values coincide. Substituting 0o = 02 = 02 



, „ ,1/3 

4 cr) 



'1 



2.9 



h 



1/3 



,29 

and this critical configuration (in the limit h — > 0) takes place at 



2 k ( 4 - 71 ) 



2^ 2/3 



(4-72) 
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for arbitrary value of to > k c . At larger values of fc the system of equations (|4.65[1 - <|4.67[1 have no solutions at all. 
The simplest way to understand this, is to consider the eas. (|4.66[) and (|4.67|) for k and to much bigger than k c . In 
this case en. (|4.fi7jl would require the values of (j>\ and <f>2 as the functions of k and to to be of order fc -1 / 2 and to -1 / 2 , 
tending to zero as m, k — > oo, while ca. (|4.()t)[) tells that both <p\ and (f>2 must remain finite. Thus, the summations in 
the instanton free energy contribution, ea. i|3.30[l . has to be limited by the finite value k c {h): 



F G (h) = -V 



(-iy 



k=a 



a-E(k,m\h) x 

. e ( ' ' _ v y 
/ .a/2 

r d (dctr) ™=fe c 

V / n=0 



(-1)' 



£3 



L^, k / „ N l/2 

k=o R D (detT) 

V / n=0 



(4.73) 



(here too = [<?/u] + 1)- This expression has to be compared with the corresponding summation at h = 0, studied in 
the previous subsection. Although all the terms in these series are non-analytic functions of the disorder parameter 
u (see previous subsection), until k,m <C k c (h), the difference between the terms with h ^ and the corresponding 
terms with h = can be represented in the form of corrections in powers of h. On the other hand, as we will 
see below, the terms with k ~ k cr (h) are non-analytic functions of h, and their differences with the corresponding 
zero-field contributions (in particular the differences of the instanton energies) can not be expended in powers of h. 
It is these terms, which are the contributions of the critical instantons which yield the non-analytic in h part of the 
free energy SFcih). Rigorous extracting of this piece must be a very difficult problem (it would require derivation of 
the instantons energy for arbitrary fc, m and h, as well as summations of the full series in ea. (|4.73(l ). but with the 
exponential accuracy the form of this non-analytic singularity is defined only by the energy E^ cr ^ (h) of the critical 
instanton, 



SF G (h) 



exp 



-E {cr \h) 



(4.74) 



Let us study the instanton configuration and estimate its energy in the vicinity of its critical state. For that let us 
rescale the fields, 



1/3 



tpi (x/R c ) 



(i = 1,2,0) where 



(4.75) 



(4.76) 



defines the spacial size of the critical instanton. In terms of the rescaled fields ip, the energy of the instanton takes 
the form 



E{k, m) 



gl/* 



d D z 



\ (Wi) 2 + ^y-^ (V^ 2 ) 2 - y (V^o) 2 + U V*. Vo) 



(4.77) 



where 



U (^1,-02,^0) 



— f [fcV'i + (m — k)ip2 — m-0 2 ] — — u [ktpi + (m — k~)^j\ — mijj^\ 



+ - [hp* + [m — fc)V>2 — m ^o] — [^"01 + (m— k)ip2 — mip ] 



(4.78) 



is the potential controlled by two rescaled parameters 



u 



gl/3 h 2/3 ' 

The instanton configuration is defined by the three equations, 

—Atpi + ftpi + — utpiS 
(i = 1,2,0) where 

S = k%jj\ + (to — k)ip2 ~ mipQ 



1 



(4.79) 

(4.80) 
(4.81) 
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Let us consider these equations in the limit h — > and and for the values of k and m of order of k c ~ h~ 2 ^ 3 — > oo. 
There are two types of terms in eas. H4.8U|) : (1) fnpi and nip^S which are both of order /i~ 2 / 3 — * oo; and (2), the rest 
of the terms, which are of order of one. In the limit h — > the diverging terms must balance each other: 

fipi ~ u (ktfjf + (m - fc)V| - rnipf) ipi (4.82) 

This condition is consistent with the requirement that both k and m are of the order of k c {h), ea. l|4.72[i . Substituting 
these estimates into eas. H4.77j l. (|4.78|l . for the energy of the instanton in its critical configuration we get 

2 

E{cr) ^r D ' 3 (4.83) 

Thus, with the exponential accuracy, the field dependent non-analytic part of the free energy has the following explicit 
form 



6F G (h) ~ e~ Eicr) ~ exp ^-( C onst)-^h- D / 3 j (4.84) 

which perfectly agree with the mean- field " hand- waving" estimate, ea. <|2.16|) . It has to be noted, however, that 
in deriving the above result we have neglected the effects of the critical fluctuations, which is justified only if we 
consider the system at temperatures not too close to the critical point, t ^> r g — t 2 /( 4 ~ d \ In the close vicinity 
of T c , at r <C T g , the situation becomes much more complicated: here one would have to combine the systematic 
(renormalization-group) integration over fluctuation with the background instanton solutions. 



V. CONCLUSIONS 



In this paper the systematic approach for the off-perturbative calculations in disordered systems has been formulated 
(Section III). It has been demonstrated how the proposed scheme works in the most simple but non-trivial case of 
weakly disordered ferromagnetic Ising model away from its critical region. Here the non-analytic (as the functions of 
the disorder parameter and the external magnetic field) contributions to the free energy has been derived, eas. 14.62fl . 
(|4.84l) . and it has been demonstrated that in terms of the replica field theory such contributions appear due to 
instanton-like excitations. Of course, it is hardly possible to register the presence of these exponentially small parts 
of the free energy both in real and in numerical experiments, and in this sense the present results has mostly pure 
theoretical interest. On the other hand, thinking about the others much more complicated problems of the statistical 
mechanics of disordered systems, the investigations made in this paper look rather promising. In particular, it does 
not look completely unrealistic to try to combine the present off-perturbative scheme with the renormalization-group 
treatment of the critical fluctuations, to settle down recent suspects^ that off-perturbative degrees of freedom could 
be quite relevant in the vicinity of the critical point, so that the nature of the phase transition in random ferromagnetic 
systems may appear to be not as simple as it was thought in early days of the theory of the critical phenomena in 
disordered materials^. 
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